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ABSTRACT 

We compare the near-infrared (NIR) H band photometric and morphological properties of low-z 
(z < 0.3) 3CR radio galaxies with samples of BL Lac object and quasar host galaxies, merger remnants, 
quiescent elliptical galaxies, and brightest cluster galaxies drawn from the literature. In general the 
3CR host galaxies are consistent with luminous (~ L*) elliptical galaxies. The vast majority of FR IPs 
(~ 80%) occupy the most massive ellipticals and form a homogeneous population that is comparable 
to the population of radio-loud quasar (RLQ) host galaxies in the literature. However, a significant 
minority (~ 20%) of the 3CR FR IPs appears under-luminous with respect to quasar host galaxies. 
All FR II objects in this faint tail are either unusually red, or appear to be the brightest objects 
within a group. We discuss the apparent differences between the radio galaxy and RLQ host galaxy 
populations. RLQs appear to require > 10 11 Mq host galaxies (and ~ 10 9 Mq black holes), whereas 
radio galaxies and RQQs can exist in galaxies down to ~ 3 x 10 10 M©. This may be due to biases in 
the measured quasar host galaxy luminosities or populations studied, or due to a genuine difference 
in host galaxy. If due to a genuine difference, it would support the idea that radio and optical active 
galactic nucleii are two separate populations with a significant overlap. 

Subject headings: galaxies: active — galaxies: fundamental parameters — (galaxies:) quasars: general 



1. INTRODUCTION 

Much effort has been invested over the last decade and 
a half to characterize the host galaxies of quasars us- 
ing the Hubble Space Telescope (HST). Unfortunately by 
their very nature quasars remain difficult to study, and 
results remain ambiguous due to selection effects. De- 
tailed isophotal analysis of quasar host galaxies is im- 
possible due to the strongly anisotropic contamination 
by the nucleus and point-spread function (PSF) of the 
instrument, which also hinder accurate spectroscopy and 
thus measurement of host galaxy dynamics. In the ab- 
sence of any forthcoming space-based coronographic mis- 
sion (e.g., the Terrestrial Planet Finder Coronograph - 
TPF-C 10 ), researchers have begun to explore the na- 
ture of type 2 (heavily obscured) active galactic nucleii 
(AG N) to study their env ironments in greater detail 
(e.g., IZakamska et all 120061 ) . Type 2 AGNs have been 
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identified fro m the Sloan Digi tal Sky Survey (SDSS) out 
to z = 0.83 (|Zakamska et al.ll2003l ). but the sample re- 
mains incomplete and biases difficult to measure. Radio 
galaxies offer us an alternative sample of type 2 radio- 
loud AGN, detectable to far higher redshifts in a statis- 
tically complete manner. 

In this paper, we compare the overall near-infrared 
(NIR; H band) galaxy properties of the z < 0.3 3CR 
radio galaxies with those of quasars and BL Lac objects 
at similar redshift, and to quiescent early-type galaxies 
and mergers available in the literature. For the reasons 
outlined above, we intend to use the radio galaxies as a 
link or bridge between the quasars and quiescent galaxies 
which can be studied in much greater structural, dynam- 
ical and population detail. Nearby radio ga laxies have 
been statistically studied in recent papers (jBest et al.l 
2005; Mau ch and Sadlenl2007t ) and are well-known to oc- 
cupy luminous galaxies in general. Indeed it is by now ac- 
cepted that the primary requirement for radio-loud AGN 
activity is a supermassive black hole > 10 9 Mq, which 
generally requires a very luminous elliptical host galaxy. 
The principal aim of this paper is to explore in detail the 
subset of galaxies drawn from the nearby universe that 
are capable of hosting a radio-loud AGN. We wish to de- 
termine what additional effects the host galaxy plays in 
determining the radio loudness of a source, and to iden- 
tify any sources that are hosted by "unusual" galaxies. 
A secondary aim is to determine whether there are any 
systematic biases introduced in the study of quasar host 
galaxies. 

Our approach is to analyze the results of isophotal 
and parametric modeling of the host galaxies of 3CR ra- 
dio so urces from | Madrid et all (|2006l) - hereafter Paper 
I, and iFlovd et al.l ( 2008f ) - hereafter Paper II, and to 



compare the results to those for similarly studied sam- 
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pies of galaxies and AGN. In Paper I we presented the 
imaging and photometry for the first part of our survey. 
In Paper II we presented the modeling used here, and 
showed that two different techniques (widely adopted by 
the galaxy morphology and quasar host galaxy communi- 
ties, respectively) reliably recovered the properties of the 
NIR host galaxies of the low-redshift (z < 0.3) 3CR ra- 
dio sources without strong optical nuclei. Furthermore, 
for synthetic quasars (created by taking a true galaxy 
and adding an artificial central point source) , the quasar 
modeling technique was shown to recover the correct un- 
derlying morphological parameters to within 10%, and 
galaxy luminosity within to 2%. 

The galaxies associated with the 3CR radio sources 
have long been known to be typically ellipti- 
cal (pvfatthcws et al.l 1 196 45. with a number of them 
associa ted with the brightest galaxi es in a cluster or 
group (|Burnsl 119901: iBest et alM2007l) . Recent surveys 
with HST have revealed a wealth of information on the 
environments of these uniquely powerful sources at a 
range of wavelengths (|de Koff et al.l 119961 : iMartel et al 
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20061 iPrivon et all 120081; iFlovd et al.l 120081 ). and uncov- 
ered n umerous new jets, dust disks, etc. Bu t until re- 
cently (jDonzelli et all 120071: Iflovd et al.ll2008D . the 3CR 
lacked a detailed and systematic classification of their 
host galaxies in such a way that can be straightforwardly 
compared with existing samples of quiescent galaxies in 
the literature. Observing in the infrared offers impor- 
tant advantages in the study of both galaxies and AGNs 
providing a direct measure of dynamically-dominant old 
stel lar population, an almost-constant mass-to-light ra- 
tio ()Zibetti et a l. 2002), and a peak in the ratio of galaxy 
to AGN continuum, giving the best chance of detecting 
the host galaxy in cases where AGN is unobscured. Thus, 
it seems a natural application of the NICMOS H band 
data set to make a direct comparison to similar sam- 
ples in the literature. We need to be certain which radio 
galaxies correspond to which RLQs in order to be able to 
use radio galaxies as a proxy for quasars in detailed spec- 
troscopic, dynamical and structural studies of the effect 
of AGN on environment and vice versa. 

This paper is structured as follows. In Section[5]we de- 
scribe the samples used in this paper, and refer briefly to 
the observations and data reduction used. In Section [3] 
we describe the results of sample-wide comparisons ac- 
cording to different observational properties of the galax- 
ies. In Section@]we discuss those comparisons, property- 
by-property across all the samples studied, and look at 
the R e — fj, e Kormendy relation, host-to-nuclear lumi- 
nosity distribution and host luminosity versus extended 
radio power distributions. We pick out several anoma- 
lously faint 3CR host galaxies for discussion in further 
detail in Section [3J In Section |H] we attempt to draw this 
discussion together and look at the issues remaining in 
the study of quasar host galaxies. Finally we conclude in 
Section[7]with a summary of our findings and suggestions 
for future study. 

2. SAMPLES, OBSERVATIONS AND DATA REDUCTION 

In this p aper we take the revi sed 3C sample (3CR) as 
defined by iSpinrad et all ()1985f ). restricted to z < 0.3, 
as a statistically complete sample of low-redshift radio 
galaxies, and compare it in various physical terms with 



samples of low-redshift quasars, quiescent galaxies and 
mergers that have deep imaging and structural modeling 
of the galaxies available in the literature. The 3CR cov- 
ers sources in the northern hemisphere that are brighter 
than 9 Jy at 178 MHz. H band snapshot observations 
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Fig. 1. — Rest-frame 178 MHz radio luminosity of the present 
sample (down-pointing triangles = FR Fs; up-pointing triangles 
= FR IPs; diamonds = HEGs; squares = LEGs), along with 
the quasars (st ars) and BL Lac o bjects (circles) with radio de - 
tections from lUrrv et all pQOOh : IMcLeod and McLeodl pOOll ): 
IDunlop et all (120031 ); IFlovd et all [120041 ). The dotted line shows 
the radio-quiet cutoff traditionally adopted in quasar studies 
(log 1 g(Li7g) = 25.83). The 3CR definition is of sources brighter 
than 9 Jy at 178 MHz. 



of th e low-z 3CR sample (as defined by ISpinrad et all 
fl985h began in HST Cycle 13, using NICMOS 2, and 
we presented the first year's imaging in paper I. Obser- 
vations continued through Cycle 14 yielding a total of 
101 sources (an 87% completeness), including 11 drawn 
from the NICMOS 2 archive. Modeling and classifi- 
cation of all 101 galaxies were presented in Paper II. 
In the present paper we revisit those modeling results, 
drawing comparisons with other galaxy samples drawn 
from the literature. See Paper I for full details of the 
3CR NIC2 observations and notes on half the sample. 
See Paper II for details of the 3CR NIC2 data reduc- 
tion, modeling results and notes on the remainder of the 
sample. In particular, the machine-readable Tables 3, 
4 and 5 of Paper II contain the basic radio properties, 
one-dimensional and two-dimensional Sersic profile fits, 
respectively, for t he entire present sample of 101 3CR 
sourc es. Also see iJenkins et al.l 119771 : iGiovannin i et all 
1 19881 for 5 GHz radio properties and iLillv and Longairl 
( 1984) for K band fluxes. All reduced, science ready data 
is available for public download from the STScI archive 
Web site 11 . 

Note that from the original sample modeled in Pa- 
per II, we omit 3C 258 from study in this paper because, 
as discussed in Paper II, it is found to be a higher red- 
shift {z > 1) RLQ located behind a z — 0.165 irregular 
galaxy. 

11 http:/ /archive. stsci.edu/prepds/3cr/ 
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TABLE 1 

Details of the comparison samples used in this paper. 



Paper 


Sample Type 




Telescope /Instrument 


Filter 


Sample Distance 


Selection Notes 


U+00 a 


BL Lac object 


30 


ffi>T/WFPC2 


i?. 


0.0 < z < 1.3 


Composite sample; 30 sources at z < 0.3 


MM01 6 


RQQ 


16 


ffi>T/NICMOS2 


11 


z rj 0.3 


Luminous QSO's at z < 0.3 


D+03 c 


RLQ, RQQ, RG 


33 


HST/WFPC2 


R 


z = 0.2 


Optically matched samples of RLQ and RQQ 


F+04 d 


RLQ, RQQ 


17 


HST/WFPC2 


1 


z = 0.4 


Optically luminous QSO's 


T+96 e 


RLQ, RQQ, RG 


33 


UKIRT /IRC AM 


K 


z = 0.2 


We use D+03 c subsample 


RJ04- 1 * 


Merger 


51 


U. Hawaii 2.2 m/QUIRC 


K 


z < 0.045 


Single-nucleus candidate disk-merger remnants 


BBF929-' 1 


Elliptical 


80 


Calar-Alto 1.23 m 


V, R, I 


D < 140 Mpc 


Nearby B band luminous ellipticals 


P99 4 


Elliptical 


3 11 


Palomar 60 inch; LCO 1 m, 2.5 m 


K 


z < 0.03 


Various nearby galaxies, groups and clusters 


M+99^' 


Elliptical 


48 


UKIRT/IRCAM3 


K 


D = 99 Mpc 


Coma cluster ellipticals 


F+06 fc 


Early type 


100 


ffST/ACS-WFC 


g,z 


D = 18 Mpc 


Virgo cluster early types 


G+96 l ' m 


BCG 


119 


CTIO 1.5m; KPNO 2.1m, 4m 


Ft 


0.01 < z < 0.054 


Abell brightest cluster galaxies 


I, 


RG 


101 


ffi>T/NICMOS 


H 


z < 0.3 


Present study 
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Note. — 1 Sample size used 

2.1. Cosmology, Photometry, and Evolution 

We throughout assume a flat, A-dominated cosmology 
in which fl\ = 0.7, and ho — 0.7. For each source we 
have fc-corrected the apparent magnitudes to rest-frame 
H band assuming a spectral index of a = 1.5 for the 
host galaxy and a = 0.2 for the nucleus where neces- 
sary (f u oc v~ a ). Gal actic extinctio n corre ctions have 
been applied following ISchlegel et al.l ()1998l ) . Compari- 
son samples have been converted to the rest-frame wave- 
length (and correct angular diameter distance) of our 
study following the same cosmological assumptions and 
using a throughput template for the NICMOS2 F160W 
(H band) filter and the appropriate filter/instrument 
configuration for the original observations. 

Note that few samples exist to precisely match our red- 
shift distribution. Out to z = 0.3 we are probing a to- 
tal look-back time of 3.4 Gyr. While evolution is not 
negligible across such a time frame, on the basis of our 
hypothesis that the 3CR represent a subsample of the 
most massive elliptical galaxies at low redshift, it is safe 
to assume that most of the mass has formed much ear- 
lier and that now the galaxies are evolving more-or-less 
passively. The radio sources themselves will have evolved 
significantly, but again, following our null hypothesis, all 
early-type galaxies should have hosted a radio source at 
some point over this epoch. For quasars we can con- 
struct samples that span a similar redshift range, while 
for quiescent galaxies existing morphological studies tend 
to concentrate on local samples, but it is informative to 
compare with quiescent galaxy samples where they are 
available. 

2.2. Comparison Samples 

We explored the literature for suitable samples 
with comparable published data morphological 
information from one- or two-dimensional fitting. 
We take BL Lac objects and quasars with deep 
imaging and well-studied host galaxies fro m the 
samples of iTavlor et all (119961,; lUrrv et all (|2000[) ; 
McLeod and McLeodl (|2001h ; IPunlop et all (|2003); 
Flovd et~ai"1 (120041, . Quiesc e nt ea rly - type g alaxie s 



Dunlop et al.| 12003) ; ^Floyd et~aT| |2004| ); iTaylor et all 119961) ; ^Rothber~ 
all U999I ); ' iFerrarese et al.l 120061) ; ^Graham et all 119961) ; " iPoetman and Lai 



nd Josepfl 130041) ; ^Bender et alj 
=r (1995); 17 Madrid ot al. 12006); 



are taken from iBender et al.l (119921,; iPahrel (|1999h ; 
llvlobasher et~ai1 (|1999f ); iFerrarese et alj {2006). Merger 



remnants we take from iRothberg and Joseph! (2004). 
Finally, we use t he brightest c l uster galaxy (BCG) 
sample studied by iGraham et al.l ill !)!)(>! . The samples 
in each of these papers are briefly discussed below. See 
Table Q] for basic details of the comparison samples at a 
glance. 

In each case, we have taken care to compare published 
properties with equivalent properties within our sample. 
A wide range of conventions are used in the literature 
for describing scale length, surface brightness and disk- 
in ess. We have made use of the conve rsions discussed 
by lMilvang- Jensen and J orgcnsenl (|1999h . In cases where 
the total magnitude of the (modeled) host galaxy was 
unavailable, we have calculated it from the basic model 
parameters using the prescriptions in the same paper. 

Note that, at present, there are insufficient quasars 
and BL Lac objects with observed host galaxies to pro- 
duce a complete or volume-limited sample. The radio 
luminosity-redshift distributions are compared in Fig- 
ure [T] The 178 MHz luminosity distribution of the com- 
bined RLQ+BL Lac object sample forms a good match 
to that of the 3CR, with a Kolmogorov-Smirnov (K-S) 
test returning D = 0.20, with p — 0.76, indicating no sig- 
nificant differences between the radio power distributions 
of the two samples. 

2.2.1. Quasars and BL Lac objects 

Urry et al. (2000): lUrrv et alj (|2000l) obtained 
HST/WFPC2 F702W {R band) images of 110 BL Lac 
objects drawn from a 132-strong sample at < z < 1.3. 
The sample is a composite of seven complete catalogs se- 
lected in the radio, X-ray, and optical. We use only those 
30 sources that fall into the redshift range of our sample 
(0.0 < z < 0.3). Urry et al. apply a one-dimensional fit 
to an azimuthally averaged radial profile of the source, si- 
multaneously fitting both host and nucleus. They adopt 
an exponential disk and a de Vaucouleurs i? 1 / 4 law, but 
do not fit Sersic models, nor ellipticities. As part of a 
survey of QSO hosts in the HST archive (D. J. E. Floyd 
et al. in preparation), we have refit the Urry et al. host 
galaxies with a Sersic mod e l, usin g the same technique 
as described in lFlovd et al.1 (|2004| ) and Paper II. We use 
these results (which are consistent in terms of host lumi- 
nosity with the original data) to compare Sersic param- 
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eters. See Figure [2] for a comparison of the Urry sample 
with the 3CR. 

Mc Leod & McLeod (2001): iMcLeod and McLeodl 
(2001) studied 16zw 0.3 optically-luminous radio-quiet 
quasar (RQQ) host galaxies in the NIR F160W filter 
(H band) using NICMOS. The authors adopted a PSF- 
subtraction technique, followed by modeling of the resid- 
ual host galaxy using both one-dimensional and two- 
dimensional analysis. We compare our two-dimensional 
results to their two-dimensional results which fit expo- 
nential disk and de Vaucouleurs elliptical models, but 
no Sersic profile. Due to the presence of saturated cores 
in almost all of the images, we were unable to remodel 
this data using our technique. Figure [3] shows the vari- 
ous properties of the sample compared with those of the 
3CR. 

Dunlop et al. (2003), Taylor et al. (1996): 
iDunlop et al.l (|2003l) studied the largest and most strin- 
gently selected sample of low-redshift quasars: 10 RLQs, 
13 RQQs, and 10 radio galaxies at z » 0.2. The radio- 
loud and RQQ subsamples were optically matched, and 
the RLQs and radio galaxies were matched in terms of 
their radio luminosities. The sample selection avoided 
core-dominated radio sources in order to minimize the 
chances of beaming producing artificially boosted ap- 
parent radio luminosities. Targets were observed using 
77ST/WFPC2 with the F675W (R band) filter. Figure g] 
illustrates the various properties of the sample compared 
with those of the 3CR. The same objects were previ- 
ously observe d from the ground using UKIRT/IRCAM 
in K band bv lTavlor etall (1 19961) - see Figure p 

Floyd et al. (2004): iFlovd et all (|200J) studied 
the host galaxies of 17 quasars at z » 0.4 (including 
those of the seven optically brightest QSO's at z < 0.4) 
using HST / WFPC2 F814W (I band) imaging. The 
lower luminosity subsample (10 objects) was divided into 
equal-sized subsamples of radio-loud and RQQs that are 
optically matched both with each other and with the 
sample of Dunlop et al. (2003). Similar care was taken 
to avoid beamed radio sources as in the Dunlop sam- 
ple. Figure [5] shows the various properties of the sample 
compared with those of the 3CR. 

2.2.2. Mergers 

Rot hberg & Joseph (2004): Rothb erg and Joseph! 
(2004) studied a sample of 51 nearby (z < 0.045) candi- 
date disk-merger remnants, all with a single nucleus, and 
based on optical morphologies to include objects with 
tidal loops, tails and shells. The study used ground- 
based K band data and one-dimensional radial profile 
fitting to elliptical isophotal values in the same way as 
our one-dimensional approach. Observations were taken 
using QUIRC at //10 on the University of Hawaii 2.2 m 
telescope. This is the only large sample or mergers that 
has been studied in such morphological detail to date. 
Figure [7] shows the Rothberg sample compared to the 
3CR. Figure |8] shows the Rothberg sample compared to 
only those 3CR sources that were flagged as mergers in 
Paper II. 

2.2.3. Elliptical and Early-Type Galaxies 

Bender, Burnstein & Faber (1992): 
iBender et al.l (|1992t ) studied a sample of local 
(D < 140 Mpc) "dynamically hot galaxies", that 



are in practice a range of spheroids and bulges selected 
in terms of their B band luminosity. The sample includes 
48 giant ellipticals (defined as having Mb < —20.5), 
20 intermediate ellipticals (-20.5 < M B < -18.5), 12 
bright dwarf ellipticals (Mb > —18.5) as well as 4 com- 
pact ellipticals, 19 bulges and 5 lower luminosity dwarf 
sp heroids. Observations and mod eling were presented 
by IBender and Moellenhofll (|1987f ) . Observations were 
made using the Calar-Alto 1.23 m telescope and a NEC 
P8603/A detector, modeling followed a similar method 
to our one-dimensional fits to elliptical isophotal values, 
presented in Paper II. The sample is compared with 
the 3CR in Figure GO We adopt this sample as a useful 
cross-section of the elli ptica l popu lation. 

Pahre (1999): iPahrd (| 19991) observed 341 nearby 
(z < 0.03) early-type galaxies in K band using the Palo- 
mar 60 inch and the Las Campanas 1 m (Swope) and 
2.5 m (Du Pont) telescopes and Rockwell NICMOS3 
HgCdTe detectors. Most are from 13 rich clusters (85%), 
but also with 12% from loose groups and 3% from the 
general field. These span the full range of early types, 
forming a useful reference for the 3CR. Modeling was 
performed in a similar way to our one-dimensional fits to 
elliptical isophotes, as described in Paper II, although all 
masking was performed automatically, instead of some 
removal by hand of obvious companion sources as in our 
modeling. No Sersic fit is given. Figure ITOl illustrates the 
various properties of the sample compared with those of 
the 3CR. 

Mobasher et al. (1999) iMobasher et al.l (|1999t ) 
studied the Fundamental Plane in a sample of 48 el- 
liptical galaxies (no lenticulars) from the Coma cluster 
(D = 99 Mpc) in K band. They were selected to have 
optical half-light diameters of less than 1 arcmin in or- 
der to fit into the field of view of IRCAM3 on UKIRT. 
Curves of growth were produced for each target, rather 
than a functional form being fitted. Effective diameter 
and effective mean surface brightness were calculated di- 
rectly from the curve of growth. Figure [IT] illustrates the 
various properties of the sample compared with those of 
the 3CR. 

Ferrarese et al. (2006): As part of the Ad- 
van ced Camera for Surve ys (ACS) Virgo Cluster Sur- 
vey, [^rrareseleF^n (|2006l ) studied a sample of 100 early- 
type members of the Virgo cluster (D = 18 Mpc) from 
dwarfs (Mb = —15.1 mag) to giants (Mb = —21.8 mag) 
imaged in F475W and F850LP (g and z bands). This 
forms a useful local early-type galaxy baseline with which 
to compare our the luminosity of our sample, allowing us 
to rank the host galaxies of the 3CR against the entire 
contents of a nearby cluster. Modeling performed by Fer- 
rarese et al. was identical to our ellipse fitting followed 
by radial profile fits to the elliptical isophote values in 
Paper II, although F errarese et al. also fit core-Sersic 
([Graham et alJ l2003h models to obje cts that are suffi- 
ciently close to detect a "Nuker"-like (|Lauer et al.lll995l) 
turnover or "core" in some of the galaxy centers. We 
compare results to their z band data, converted to H as 
described above. Figure [T2l illustrates the various prop- 
erties of the sample compared with those of the 3CR. We 
note that many of these objects have noticeable disks in 
addition to their bulges (SO types), reflected in their in- 
termediate Sersic index values. This use of a single Sersic 
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Fig. 2. — Histograms for the present 3CR sample (open) and the Urrv ct al. (2000) BL Lac object sample (shaded) for the following 
properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity; Sersic index. L* is indicated by a 
dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 3. — Histograms for the present 3CR sample (open) and the McLcod and McLcod (2001) RQQ sample (shaded) for the following 
properties (from left to right): host galaxy luminosity (converted to H band); scale length; ellipticity. L* is indicated by a dotted line on 
the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 4. — Histograms for the present 3CR sample (open) and the Dunlop ct al. (2003) quasar sample (full sample - shaded; RLQs - 
black) for the following properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity; Sersic index. 
L* is indicated by a dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 

profile in such cases is questionable. 

2.2.4. BCGs 

Graham et al. (1996): iGraham etaLI (|1996l) took 
the sample of 119 Abell BCGs (0.0 1 < z < 0.054) origi- 
nally o bserved in Cousins R band bv lPostman and Lauerl 
(fl995T) on the CTIO 1.5 m and the KPNO 2.1 m and 4 m 
telescopes. Sersic profiles we re fit to the semi-ma j or axi s 
surface brightness profiles of iPostman and Lauerl (|1995h . 
similar to our one-dimensional fitting of profiles to the el- 
liptical isophotal levels in Paper II. Figure H3] illustrates 
the various properties of the sample compared with those 
of the 3CR. This remains a state-of-the-art data set for 
the study of BCG's in spite of the age of the observations. 

3. RESULTS 

For the full tabulated results of the morphological mod- 
eling of each individual source, see Paper II, Tables 4 and 
5. Here we only discuss those results in a sample- wide 
comparison. Histograms comparing the distributions of 
various properties of these samples with those of our 3CR 



sample are presented in Figures I2HT31 We performed a 
K-S test to compare each property of each comparison 
sample to the equivalent property of our 3CR sample. 
This returns the K-S statistic, D, reflecting the cumula- 
tive difference between two samples, and a likelihood, p 
that the two samples are drawn from the same popula- 
tion. Results of the K-S test, showing the goodness of 
fit between each sample and the 3CR are presented in 
Table Table |3] shows each sample compared to just 
the FR I sources in our 3CR sample. Table 0] shows each 
sample compared to just the FR II sources in our 3CR 
sample. In the following we present a broad sample-by- 
sample comparison, followed by more detailed compar- 
isons in each observational metric (morphology, elliptic- 
ity, luminosity, scale length and the R e — fi e Kormendy 
relation). In the discussion Section 2] below, we engage 
in a more detailed discussion of each property. 

3.1. Sample-by-sample Comparisons 

The BL Lac object host galaxies of the lUrrv et all 
(2000) sample span the upper luminosity and scale length 
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Fig. 5. — Histograms for the present 3CR sample (open) and the lTavlor et al] (1996) quasar sample (full sample - shaded; RLQs - black) 
for the following properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity. L* is indicated by a 
dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] The Taylor et al. 
sample is the same as the Dunlop et al. one, but studied at K band from the ground (see the text). 
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Fig. 6. — Histograms for the present 3CR sample (open) and the lFloyd et al . (2004) quasar sample (full sample - shaded; RLQs - black) 
for the following properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity; Sersic index. L* is 
indicated by a dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 7. — Histograms for the present 3CR sample (open) and the Rothbcrg and Joseph (2004) merger sample (shaded) for the following 
properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity; Sersic index. L* is indicated by a 
dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 

range of the 3CR host galaxies - Figure [2j Galaxies as 
dim (small) as the very bottom end of the 3CR galaxy 
luminosity (scale length) distribution are conspicuously 
missing, and this is also found with the other active 
galaxy samples discussed below. We return to discuss 
this point in some detail in Sections 14.21 and [SJ Interest- 
ingly, the BL Lac object ellipticities overlap nicely with 
the round majority of the 3CR, but lack the eccentric tail 
representing merging sources among the radio galaxies. 
Overall the BL Lac objects represent a poor fit to the 
3CR sample. Even considering the FR I and FR II sub- 
samples separately produces little improvement in the 

match; 

The lMcLeod and McLeodl (|2001ft RQQ host galaxy lu- 
minosities and ellipticities exhibit a good match with the 
3CR - Figure G3 The quasar host galaxies appear to be 
somewhat smaller on average, in contrast to the other 
samples of quasars discussed below that include both 
radio-loud and radio-quiet sources. Redshifts range from 
0.15 to 0.4 with an average at z = 0.25. 



The iDunlop et al.l (|2003t ) quasar host galaxies appear 
to be consistent with the bright end of the 3CR host 
galaxy distribution, assuming a simple conversion from 
RtoH b and - see Figure B} The more certain conversion 
from K (jTavlor et al.|[l9 96) to H provides a similar pic- 
ture (Figure [5]) though based on less certain host galaxy 
luminosities due to the less stable ground-based PSF. 
The higher ellipticities and scale lengths in the Dunlop / 
Taylor quasars are again possibly due to PSF correction 
problems. Note that the Dunlop et al. sample exhibits a 
strong cutoff in Sersic index in their radio-loud subsam- 
ple that is not mirrored in the 3CR population. 

The quasar host galaxy sample of iFlovd et al.l (|2004[ ) 
again matches the upper reaches of the 3CR host galaxy 
luminosity distribution, while the scale length, ellipticity 
and Sersic indices all show a similar spread in values - 
Figure [5] Note that a key finding of both Floyd et al. 
and Dunlop et al. was the relative similarity of RLQ and 
RQQ properties, although the smallest and dimmest of 
the quasar host galaxies in each sample are radio quiet in 
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Fig. 8. — As for Figure [7] but for Rothberg & Joseph (2004 - shaded) compared to only those 3CR sources that are : 
remnants (open). 
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Fig. 9. — Histograms for the present 3CR sample (open) and Bender ct al. (1992) early-type galaxy sample (shaded) for the following 
properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity. L* is indicated by a dotted line on the 
left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 10. — Histograms for the present 3CR sample (open) and the Pahrc (1999) early-type galaxy sample (shaded) for the following 
properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity. L* is indicated by a dotted line on the 
left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 



each case. This is reflected in the histograms presented 
in Figures @] and [5] 

In general, and unsurprisingly, the Rothberg & Joseph 
(2004) sample forms a poor match to the 3CR host galax- 
ies (Figure [7]). However, if we examine only the sub- 
set of 3CR host galaxies that are classed as exhibiting 
merger signatures (as defined in Paper II), we see some 
improvement (Figure [5]). In particular, there is a notice- 
able increase in the likelihood that the two samples are 
drawn from the same ellipticity distribution. But even 
among these merging 3CR sources, the host galaxies are 
significantly less disky than the general merger popula- 
tion explored by Rothberg & Joseph and the disk-merger 
population contains many sources much smaller than the 
3CR mergers, while the luminosity distributions span a 
simila r range. 

The iBender et al.l (|1992f ) sample forms a poor match 
to the 3CR host galaxies in all observational metrics (see 
Figure [9]), being generally larger and more luminous, as 
well as slightly more eccentr ic in shape. 

Generally the lPahrel fl999) sample forms a good match 



to the 3CR host galaxies, although the Pahre sample does 
exhibit a long low-luminosity / compact galaxy wing, 
beyond the cutoff in the 3CR host luminosity and scale 
length distributions - Figure QUI In terms of ellipticity, 
the m atch is excellent. 

The [Mobash er et al.l (|1999t ) sample forms a very good 
match in terms of host galaxy luminosity to the 3CR - 
Figure [TT] It is noticeable that the 3CR galaxies tend to 
be slightly larger, with a poorer match between the scale 
length distributions, due to the field-of-view imposed size 
cutoff in the Mobasher sample. 

Unsurprisingly, the 3CR host galaxies are much larger 
and more lum inous than the typical cluster galaxy repre- 
sented by the iFerrarese et al.l (|2006t ) sample - Figure [T2l 
However, they also exhibit larger Sersic indices in gen- 
eral, in spite of the fact that Ferrarese et al. deliberately 
study only morphologically-selected early-type galaxies. 
Many of the Ferrarese objects are SO's with noticeable 
disks in addition to their bulges, yielding intermediate 
Sersic values. SO's are not accurately represented with a 
single Sersic profile. 
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Fig. 11. — Histograms for the present 3CR sample (open) and the Mobasher ct al. (1999) early-type galaxy sample (shaded) for the 
following properties (from left to right): host galaxy luminosity (converted H band); scale length. L* is indicated by a dotted line on the 
left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 12. — Histograms for the present 3CR sample (open) and the Fcrraresc ct al. (2006) galaxy Virgo cluster galaxy sample (shaded) 
for the following properties (from left to right): host galaxy luminosity (converted H band); scale length; ellipticity; Sersic index. L* is 
indicated by a dotted line on the left-hand figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 
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Fig. 13. — Histograms for the present 3CR sample (open) and the lGraham et al. ( 1996) BCG sample (shaded) for the following properties 
(from left to right): host galaxy luminosity (converted H band); scale length; Sersic index. L* is indicated by a dotted line on the left-hand 
figure. The results of a two-sided K-S test for each distribution are presented in Table [2] 



The 3CR host galaxies sp an a far lower r ange of lu- 
minosity than the BCG's of Graha m et al.1 (|1996l ). are 
significantly smaller and have lower Sersic indices - Fig- 
ure 1131 The distributions have dissimilar shapes but 
there is significant overlap in terms of Sersic index. 

3.2. Host Galaxy Morphology 

We find a broad distribution of Sersic indices in our 
sample, spanning the full range from n — 1 (exponential 
disks) to n = 4 (de Vaucouleurs ellipticals), and higher. 
In Paper II we discussed the FR II radio sources that are 
hosted by disk-like (n < 2) galaxies, all of which have 
close companion objects and are either merging, or are 
likely to be post-merger systems. In comparison with 
normal cluster early- types (Figures [91 -H"2"j). this range of 
Sersic index is not unu sual. For examp l e, by comparison 
to the local sample of iFerrarese et all (|2006l ). the range 
of Sersic index is identical, but the peak in Sersic index 
in our sample falls significantly higher: in the Virgo clus- 
ter the peak Sersic index is at n — 2, with a significant 
lower- n tail, whereas in the 3CR, the peak is at n = 3, 



and the tail exhibits a much more rapid falloff to low- 
n (Figure [T2"|) . When we compare the Ferrarese sample 
to only the FR I's in the 3CR there is a dramatic im- 
provement (D — 0.10, p — 0.99), with both populations 
exhibiting similar peaks in the Sersic index distribution. 
However, whereas the Ferrarese objects are typically SO's 
with distinct bulge and disk components, the FR I's are 
typ ically found in diffuse giant (D, cD) elliptical galax- 
ies ([Owen and Laing 1989). Both populations are better 
described using two Sersic profiles. 

In comparison with the BCGs of iGraham et~aT1 (fl996h 
there is again a clear distinction, with the BCG's exhibit- 
ing a far larger range of n to high values (Figure [T5|) . 
The peak in the BCG Sersic index distribution is close 
to n — 5, with a number of objects having n > 10. How- 
ever, at the low-n end, the 3CR closely follows the cut- 
off in Sersic index seen in these most massive cluster- 
dominating galaxies. These source s are also be better fit 
by a two-Sersic component model (|Seigar et al. 2007). 

Morphology also distinguishes the 3CR host galaxies 
from typical merger systems that have Sersic indices even 
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lower than the cluster ellipticals, peaking at n = 1 (Fig- 
ure [7J) . These sources are posited as post-merger sys- 
tems on their way to forming massive ellipticals, but 
not g enerally the first-ranked ell i ptical galaxies in clus- 
ters ([Rothberg and Joseph! 120041 : iRothberg and Fischerl 
12010( 1 . Even when comparing to only the 3CR sources 
nagged as mergers in Paper II, we see a far stronger 
tendency for the 3CR mergers to be identified as bulge 
dominated (n > 2) than the general merger population 
(Figure©. 

The morphologies of the sampl e do fit in we l l with 
the (far smaller ) QSO samples of iDunlop et al.1 (2003) 
and lFlovd et~all (|2004) - see Figures SHU However, these 
QSO hosts appear to show a bimodal distribution, with a 
smaller disky (n = 1) peak, as well as the ellip tical peak 
close to n = 4. As noted bv lFlovd et al.l J2004), however, 
the three disky objects are all radio-quiet quasars, at low- 
optical luminosity similar to a classical Seyfert galaxy. 
The same is true for the Dunlop sample. Thus, while the 
3CR are morphologically similar to the RLQs studied in 
each of these papers, they are distinct from the general 
host galax y populatio n of R QQs. Using our own Sersic 
fits to the lUrrv et al.l (2000) BL Lac object sample we 
find a similar spread in Sersic indices, although there is 
a long tail to high values among the BL Lac objects. 

3.3. Ellipticity 

The ellipticity distr ibution of the 3CR matches very 
well with that of the iPahrd (|1999t l sample of K band 
studied elliptical galaxies drawn from a range of en- 
vironments (K-S t est: D = 0.06, p = 0.99). How- 
ever t he samples of iBender et al.l (|1992l) ; iFerrarese et al.l 
(|2006f) . which focus on entire populations from a clus- 
ter, show (typically E2-E3) ellipticity distributions that 
are closer to that o f the mergers than that of the 3CR. 
iDisnev et~a l. (1984) also found that radio-loud ellipticals 
are inherently rounder than radio-quiet ones. Unfortu- 
nately, no ellipticity data is available for the BCG's, but 
they are quite round in general (from visual in spection 
of the original data of lPostman and Lauerlfl995Tl . 

The quasar data seems to support a similar ellipticity 
distribution for both the Floyd and the Dunlop samples 
{D = 0.22, p = 0.42 and D = 0.14, p = 0.67, respec- 
tively). In particular, the radio- loud objects exhibit sim- 
ilar peaks at around El in both samples, but there is 
an additional highly-elliptical (E5) peak in both quasar 
samples that may be due to small number statistics, or 
indicate objects in which asymmetries from the PSF have 
interfered with proper characterization of the quasar host 
galaxy. 

The most striking cont rast in terms of elliptici t y is w ith 
the merger remnants of IRothberg and Joseph! (|2004f l - 
see Figure [7] (D = 0.34, p = 0.001). The merger 
remnants typically have much higher median ellipticities 
(E3) than those of the present 3CR sample (El). This is 
perhaps unsurprising given that the Rothberg & Joseph 
sample is targeted at objects that have recently under- 
gone a major merger and have not had time to fully relax 
to their (presumably elliptical galaxy) endpoints. But 
it does demonstrate quite graphically that the powerful 
radio sources are not, in general, consistent with recent 
major-merger systems. However, if we turn our attention 
to those objects in the 3CR that exhibit any kind of tidal 
or merger-like disturbances (flagged with mergers, shells 



or tidal tails in Paper II), we find a somewhat different 
situation - see Figure|8] a bimodal ellipticity distribution 
that follows the ellipticity distribution of Rothberg & 
Joseph's mergers somewhat better (D = 0.20, p = 0.22). 

3.4. Galaxy Luminosity 

The low-z 3CR host galaxy population peaks at close 
to L* (Mh = —23.9), but with significant spread. The 
high end is dominated by the FR II sources, but there 
are several FR II sources in the low-luminosity wing as 
well. The merger systems studied by Rothberg & Joseph 
are generally as luminous as the 3CR host galaxies, if not 
more so. Rothberg & Joseph conclude that these systems 
end up largely as L* (or brighter) elliptical galaxies after 
the starburst fades. 

Comparing with elliptical galaxies in general, we see 
a good m atch in host luminosity distribution with the 
sample of lMobasher et al.l (|1999t l (D = 0.11; p = 0.82) 
- particul arly to the FR I's {D = 0.12; p = 0.83). The 
sample of lPahrd (|1999h gives a somewhat poorer match 
(D = 0.15; p — 0.11) which is improved if we look only 
at the FR IPs {D = 0.15; p = 0.26). The IFerrarese et all 
(2006) sample covers the entire Virgo cluster, and un- 
derstandably peaks at far lower luminosities, but it is 
instructive to see the crossover between the 3CR and the 
entire contents of such a cluster (Figure IT2|) . The 3CR 
hosts are representative of the top one third, in terms of 
luminosity, of the galaxies in Virgo. While there is some 
overlap with the BCG population, in general when com- 
pared the 3CR are smaller and less luminous. While the 
FR I's were found to have a similar Sersic index distribu- 
tion, any apparent similarity with the Ferrarese sample 
is eliminated by comparing their luminosities. 

The median luminosity of the 3CR host galaxies is 
somewhat dimmer than that of the typical quasar host 
galaxy, or BL Lac obj ect, based on simp l e transforma- 
tions from the fluxes of Urrv et aTl (|2000h : IDunlop et al.1 



(|2003f) and iFlovd etal 



(2004D assuming an old stellar 
2.11) . We acknowledge that there 



population (see Section 
are difficulties in comparing to luminosities measured in 
different bands in this manner, but with the existing sam- 
ples of low-redshift quasar host galaxies, a single band 
is generally all that can be compared, and so there is 
no constraint on the mass-to-light ratios of these ob- 
jects. However, for the Dunlop et al. sample, ground- 
ba sed K band imagin g was published in an earlier paper 
bv lTavlor etafl (|1996l ). offering us a convenient test-case: 
although the K band morphologies are less reliable due 
to the problem of disentangling the host and the nucleus 
from the ground, luminosity in most cases is well con- 
strained, and they can be much more reliably converted 
into an H band luminosity. We see an essentially identi- 
cal luminosity distribution to that exhibited by the Dun- 
lop/Taylor data - Figures |4] and [5] Once again, we see 
a significant population of sub-L* 3CR hosts that are 
not mirrored in the quasar population. As was stated 
earlier and in Paper II, the low-luminosity 3CR objects 
tend to be low-re dshift, low-radio-p ower FR I sources. 
We note that the lUrrv et all ((2000) sample exhibits a 
slightly greater crossover in luminosity with the 3CR at 
this low end, in accord with the unification of these ob- 
jects with the BL Lac objects, but there remains a deficit 
of the faintest objects seen in the 3CR. 
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3.5. Galaxy Scale Length 

The 3CR host galaxy scale lengths, in keeping with 
the host galaxy luminosities, are typical of large (~ L*) 
elliptical galaxies, with a median of 3.6 kpc, and a few 
sources larger than 10 kpc in size. They are somewhat 
smaller than reported in many parts of the quasar litera- 
ture, but are consistent with those of iFlovd et al.l ([2004D 
(D = 0.19, with p = 0.65) and somewhat consistent 
with iMcLeod and McLeodl l| 20011 ) (D = 0.32, with p = 
0.09). They are inco nsistent with the s cale length distri- 
butions rer20rted_bxJUrr^eraD (|200Q[ ) (D = 0.33, p = 
0.002). IDunlop et al.l (120031) (D = 0.52, p ~ 10~ 4 ) and 
iTavlor et all (|l996D (D = 0.63, p ~ 10~ 9 ). The latter 
data set in particular is known to overestimate the scale 
lengths, due to atmospher i c seei ng problems: compare 
results with IDunlop et al.l (|2003f ) , who study the same 
sample in R band using HST data - Figures [4] and [5l The 
Taylor et al. ground-based work g enerally recovers acc u- 
rate host galaxy luminosities (see IDunlop et all ((2003)), 
but poorly constrained morphological parameters, such 
as scale length and surface brightness. That the 3CR 
scale lengths agree at least with the Floyd et al. quasar 
sample is reassuring, and it is possible that earlier studies 
overestimated the quasar host galaxy scale lengths some- 
what. Considering the 3CR FR II sources alone, we get 
a considerable improvement to the comparison with the 
Floyd et al. quasars (D = 0.18, p — 0.78). Considering 
3CR FR I sources alone we find a significant improve- 
ment in the comparison to the IMcLeod and McLeodl 
(|20Qll ) RQQ host galaxy scale lengths (D = 0.13, with 
p = 0.99). 

While the Merger sample of Rothberg & Joseph (2004) 
is generally as luminous as the 3CR, we find that in gen- 
eral, the 3CR host galaxies are somewhat larger, consis- 
tent with the merger's fluxes being dominated by active 
star formation in their inner regions. 

Comparing to normal elliptical galaxies there are no 
statistically good matches in terms of scale length. The 
3CR are slightly larger on average than the Pahre sample 
and Mobasher et al. sample of ellipticals, but smaller 
than the BBF92 sample. Interestingly, the 3CR span 
the noticeable gap in the Ferrarese sample scale length 
distribution (Figure [12]). The 3CR are smaller in general 
than the Graham et al. BCG's which lack objects smaller 
than 2 kpc. 

3.6. Kormendy Relation (R e -fi e ) 

The strong correlation between scale length and sur- 
face brightn ess is a well-know n feature of bright ellipti- 
cal galaxies (|Kormendvlll977l ). and is a projection of the 
more physically significant "fundamental plane" . The re- 
lation followed by the present low-redshift 3CR sample 
is close to that of quiescent massive elliptical galaxies, 
following a slope of 2.97±0.1 (Paper II). The BCG's fol- 
low a similar trend but extend to far larger radii. The 
same rel ation is also a g ood fi t to th e quasar hos t galax y 
data of IDunlop et ail (|2003l ) and IFlovd et al.1 (|2004D . 
This provides strong evidence that all of these objects 
are structurally similar, dynamically relaxed systems, al- 
though proof of this requires spectroscopy to study the 
dynamics. Mergers too fit remarkably well with the over- 
all trend in the AGN and quiescent elliptical popula- 
tions. Putting all of the samples together, including the 



Virgo cluster galaxies of iFerrarese et al.l (|2006l) results 
in the plot shown in Figure El It is clear here that 
the giant ellipticals, and powerful AGNs (quasars, radio 
galaxies, BL Lac objects) form a single population in the 
R e — n e space. Dwarf ellipticals do not obey the same 
trend (jCapaccioli and Caonlll99H iGraham and Guzman! 
2003)): for a given scale length these objects are far dim- 
mer than their cousins on the Kormendy relation, and 
thes e objects are clearly not rep resented by the 3CR. 
The iRothberg and Josephl (|2004l ) mergers appear to be 
on the low side of the typical elliptical galaxy Kormendy 
relation. 

4. DISCUSSION 

4.1. Host Galaxy Morphology 

The modeling used in this paper is deliberately sim- 
ple, aimed at fitting the bulk of each galaxies flux with 
the minimum of parameters. Many if not all luminous 
galaxies have core regions that deviate from simple Sersic 
models, and which may be detected in many of the clos- 
est objects presented here, but in the interests of com- 
paring like with like we have adopted a Sersic model and 
compared to the large literature on fitting similar mod- 
els to control samples of galaxies and quasars. Regard- 
ing morphology of quasar host galaxies it is difficult to 
distinguish an n = 6 bulge, say, from an n — 4 one 
underneath the glare of a quasar nucleus (My < —23). 
The spread in n present in the 3CR host galaxies does 
in fact connect them more firmly with the quiescent el- 
liptical galaxy population, which does not in general fol- 
low a perfect de Vaucouleurs type r 1 / 4 law which tends 
to be seen in quasar host galaxies. It seems likely that 
when applying such morphological models to the host 
galaxies of quasars, there is insufficient signal under the 
photon shot noise in the central region to accurately de- 
termine n in many cases: although it is clear that in 
most cases the difference between an elliptical or bulge- 
dominat ed and disky or d i sk-dominated host is mea- 
surable dMcLure et al.lll999t IMcLeod and McLeodll2001l : 
IDunlop et al.l l2003t IFlovd et al.ll200l . it is difficult to 
determine the precise morphology as characterized by a 
single continuous parameter given current observational 
technology. 

The 3CR host galaxies follow the same Kormendy rela- 
tion to the host galaxies of quasars, and giant ellipticals 
(Figure I14|) : even the very lowest luminosity ones are 
clearly distinguished from the dwarf ellipticals. However, 
what is clearly missing is dynamical information on these 
sources that would provide a real physical insight into the 
conditions inside the host galaxi e s. A study of southern 
radio galaxies by iBettoni et al.l ()200lD has shown that 
those sources are indeed hosted by dynamically relaxed 
elliptical galaxies. The closest comparison sample is that 
of Mobasher et al. (1999) which already has suitable 
spectroscopic data. It would be a valuable study to in- 
vestigate the relative dynamical states of these two sam- 
ples. Spectroscopic data for the 3CR host galaxies would 
enable us to determine the dynamics and stellar popu- 
lations of these objects in order to explore their ages, 
merger histories, and the differences between sources ly- 
ing in rich environments and those lying in the field. 

4.2. Host Galaxy Luminosity 



Radio galaxies in context 



11 



C\2 
I 



o 

CD 
O 

S-H 

cd 

CO 



14 



16 - 



18 - 



20 - 



22 - 



24 - 



26 




+ 3CR 

*£f Quasars, BL Lacs 

O Ellipticals 

□ BCGs 

A Mergers 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 



lo gio( R i/ 2 / k P c ) 



Fig. 14. — scale length-surface brightness for the 3CR along with comparison samples (after conversion to H band) of quasar and BL Lac 
object host galaxies, quiescent ellipticals, BCG's and mergers as presented in the text. Best fitting scale length-surface brightness relations 
are shown for the FR I's (dashed) and the FR IPs (dot-dash). The locus of an L* galaxy is shown by the dotted lin e. The 3CR fits i n wel l 
with the structural properties of the giant ellipticals, and is distinct from the dwarf elliptical subsample within the Fcrrarcsc ct al. (2006) 
sample from the Virgo cluster. 

The most surprising finding from this survey has been 
the relative faintness in H band of a small number of the 
3CR host galaxies relative to the luminosities expected 
from studies of quasar host galaxies. Even when we con- 
sider only FR II sources, several host galaxies exist that 
are rather fainter than the cutoff of ~ L* inferred from 
the RLQ host galaxy luminosity functions of Dunlop ct 
al. (2003) and Floyd et al. (2004). Assuming a mass-to- 
light ratio for early- type gal axies in general, Ml Lh ~ 1.2 
(with a scatter of 0.2 dex IZibetti et al.l 120021 ) we have 
just three sources with masses > 5 x 10 11 M Q : 3C 130, 
3C 196.1 and 3C 338. Firstly we note that the me- 
dian mass of the sample estimated in this simplistic 
way is 10 11 M Q and the mean mass 2 x 10 11 M . 
This is significantly higher than the ~ 3 x 10 10 Mq 
mass cutoff inferred for a transition in galaxy formation 



mode (jKauffmann et al.l 120031: iKhochfar and Silkl l2006l 
|2009| ). Secondly, the vast majority of the anomalously 
faint tail (Jkfg > —23) in the 3CR host galaxy distri- 
bution is made up of very nearby (11 at z < 0.05, 19 
at z < 0.1) radio sources that are relatively radio-dim, 
some even classically "radio-quiet" by the simplistic stan- 
dard adopted for the RLQ/RQQ break in quasar studies. 
These sources would not be selected as part of the 3CR if 
they lay just a little more distant and cannot properly be 
compared with RLQs. Just two of the sources with faint 
host galaxies, 3C 61.1 (an FR II) and 3C 314.1 (an FR I) 
lie at higher redshifts (z — 0.186 and z — 0.119, respec- 
tively). In addition, seven of the lower redshift sources 
are worth considering: 3C 105, 3C 198, 3C 227, 3C 326, 
3C 353, 3C 402 and 3C 445 are all FR II radio sources 
in abnormally dim galaxies. These objects are discussed 
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in greater detail in the following section. Finally, we 
remark that just two of these dim FR II host galaxies 
(3C 61.1 and 3C 445) have simple estimated masses be- 
low the ~ 3 x 10 10 Mq mass cu toff inferred for a tran- 
sition in galaxy formation mod e (jKauffmann et al"1l2003l : 
iKhocMar anTSiigi20M 120091) . 

It is interesting to explore how this finding compares 
to statistical surveys of nearby radio galaxies. The 
largest sample is that of iMauch and Sadlerl (|2007l ) who 
examined 7824 radio sources from the 1.4 GHz NRAO 
VLA Sky Survey (NVSS) with galaxies brighter than 
K = 12.75 mag from the 6 degree Field Galaxy Survey 
(6dFGS). The galaxies in the 6dFGS span the redshift 
range 0.003 < z < 0.3, making an ideal comparison to 
the samples presented here. Typical colors for early- type 
galaxies are H — K » 0.2 — 0.5. Galaxies with luminosi- 
ties Mk > —23 are extremely rare - rarer than we find 
for Mh > —23. However, we note that the 6dFGS cut- 
off will lose Mk > —23 galaxies at z ss 0.033 and even 
Mk > —23.5 galaxies at just z w 0.041. Thus, it appears 
that this rare population was essentially invisible to ear- 
lier surveys and have only been discovered by virtue of 
deep imaging of a purely radio-selected sample. 

4.3. FR-type and Unification 

iFanaroff and Rilevl (| 19741 ) showed that radio galax- 
ies exhibit a change in morphology at a monochromatic 
power corresponding to roughly 10 24 ' 5 W Hz -1 at a rest- 
frame frequency of 1. 4 GHz (convert ed from the original 
cutoff at 178 MHz - iBicknelll [1995D . It is understood 
that the two clas ses have differ ent bulk energy trans- 
port mechanisms (lBicknellll995l and references therein). 
lOwen and Ledlowl ( 19941) showed that on the optical- 
radio plane, the division between the two classes, which 
occurs o ver a full decade i n radio power, becomes much 
cleaner. lUrrv et ail (|2000l) took this a stage further by 
looking at just the extended radio power (subtracting off 
the core from the total) and the optical luminosity of the 
host galaxy alone. To do this, they used the 5 GHz ra- 
dio flux which is sensitive to core emission and provides a 
core/extended distinction for the largest possible sample. 
This avoids the effects of beaming and cleans up the di- 
vision between FR I / BL Lac object and FR II / quasar 
even further. The dividin g line was found to follow that 
predicted by the models of lBicknelll (|1995[ ). with the FR I 
sources containing turbulent, transonic jets, a nd FR II 
sourc es jets that are super or hypersonic (see IBicknelll 
1 19951 and references therein). This line was later found 
to be consistent with a constant Eddington accretion 
ratio ()Ghisellini and Celottil |2001| ). In Figure [15] we 
plot the available 5 GHz data for the presen t sam- 

alone 



plot tne available o <j j 1z data tor tne presen i 
pie (| Jenkins et al l 119771: iGiovannini et all 11988 h 



withRLQs from iDunlop et al.l (120031) andlFlovd et all 
(|2004l) . and BL Lac objects from lUrrv et~all (|2000l) . The 
division confirms once again that FR I's are unified with 
BL Lac objects, while FR IPs are analogous to RLQs. 

We also note that a number of "radio-quiet" quasars 
studied by Dunlop et al. and Floyd et al. may actually 
occupy a similar region of the host luminosity - extended 
radio power plane to the FR I's and BL Lac objects, if 
their radio luminosities are dominated by an extended 
component. This is uncertain since RQQs have not been 
studied in such great depth in the radio as RLQs, and 
nothing is known of their radio structure. 



5. SOURCES WITH FAINT HOST GALAXIES 

In this Section we discuss the eight FR II radio 
sources with anomalously faint (Mh > —23) host galax- 
ies and the unusually radio-powerful FR I, 3C 314.1. 
Extinction maps were generated for thes e nine sources 
by d i viding the optical WFPC2 images (|de Koff et al.l 
[19961 iMartel et alJU99flj ) by the NICMOS images, after 
registering, scaling and smoothing the former to the 
plate scale and resolution of the latter. In several cases, 
the optical-IR colors suggest strong extinction, sufficient 
to close some of the gap to L* in H band, but not all of 
it. Note there are an additional four with intermediate 
luminosities -23 < M H < L*: 3C 111, 3C 234, 3C 371 
and 3C 390.3 which we do not discuss here. 

3C 61.1 is an FR II radio source in an under-luminous 
host galaxy [Mh = —22.2), with apparent spiral arm 
structure (see paper II and references therein) in spite 
of its high radio power (log 10 ^178 = 26.27). However, 
this source is the brightest member of a small group of 
galaxies visible on the NICMOS image, and is probably 
associated with a larger clu ster that includes a radio- 
quiet AGN at z = 0.184 (|Kristian et al.l 119781) some 
30" to the ea st. It is clearly very dusty from the optical 
morphology (|de Koff et al J 1199a ). The extinction map 
is complicated by the star formation in the apparent 
spiral arms, but the central region appears to be heavily 
obscured, and we estimate ~ 0.5 H mag of extinction 
in the central reg ion of this so urce (assuming a diffuse 
dust model after iMathisI 1 19901) . giving an unextincted 
H band luminosity of Mh — —22.7. It is likely to be 
the brightest cluster member (though this is uncertain 
without a proper study of the host galaxy of the nearby 
AGN). 

3C 105 is an FR II radio source at z = 0.089, hosted 
by a highly flattened and faintly disturbed (in our H 
band image) elliptical with M#(Host) = -22.88. The 
host galaxy of 3C 105 is known to be e xceptionally red, 
with strong stellar abso rption lines ([Tadhuntcr et al.1 
Il993t ICohen et al.1 I1999D . and both existing optical 
HST images provide insufficient detection to produce a 
color-map. We estimate an R — H color for the host 
galaxy of ~ 8 ± 2 (very approximate due to the barely 3ct 
detection in R band) indicating a very heavily obscured 
host galaxy. Cohen et al. commented that the lack of 
polarized flux could be due to the obscuration of the 
scatterers, but discount this hypothesis on the basis 
that the obscuring material would have to cover a large 
fraction of the ISM. The HST images, although shallow 
in the optical, suggest that a large quantity of absorbing 
material could indeed be the cause. 
3C 198 is an FR II radio source at z = 0.081 
in a slightly elongated elliptical host galaxy with 
Mff(Host) = —22.55 and a detectable nuclear point 
source. The IR nucleus and lack of common reference 
points on both the H band and optical images mean 
that we were unable to produce a reasonable color-map. 
This source app ears to be the brightest member of a 
small group (see lMcHard^fl97l . 

3C 227 is a broad line FR II radio galaxy at z = 0.086. 
We can clearly see the host galaxy (Mh = —22.66) and 
the bright active nucleus on this NIR image. Both the 
optical and UV images show a very bright unresolved 
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Fig. 15. — FR type diagnostic can be made according to optical luminosity vs. radio power, according to IQwen" and LedlovJ $199$) . 
Using only extended radio power and host galaxy luminosity avoids the effects of beaming l|Urrv et al.N2006T ). Models from Bickncll (1995) 



are overplotted (solid line: 7max/7n 



10 4 , v c = 10 1U Hz, / = 1; dotted line: 7 max /7n 



10 4 



10 Hz, / = 0.5). For comparison 



we have plotted the RQQs (filled stars) from Dunlop et al. (2003) and Floyd et al. (2004), assuming that all their radio flux is extended. 
Symbols used are the same as in Figure [T] with filled symbols indicating radio-quiet objects and open symbols radio-loud objects. 

the optical-IR color-map, and is unusually blue for an 
elliptical gala xy (R — Hk. 1.5). We no te that this object 
is flagged by Buttiglio ne et al.l ()2009f ) as a relic radio 
galaxy which is no longer being actively powered by the 
jet. 

3C 326 is an FR II radio sourc e and exhibits a LINER 
spectrum (|Simpson et al.l I1996D . and the host galaxy 
is quite dim {M H = -22.9 9) and red (R - H = 2.4). 
From the optical image (|de Koff et al.1 fl99l . it is 
extremely dusty, and appears to be an edge-on SO. The 
color-map shows up to 0.3 H mag of extinction in the 
centr al regions of 3 C 326 (assuming a diffuse dust model 
after iMathis! |l990j) , suggesting an unextincted galaxy 
luminosity of Mh ~ —23.3. 

3C 35 3 is an FR II radio source with a LINER spec- 
trum (jSimpson et al.l [1996) hosted by a round giant 



nucleus (|Martel et al.l 119991: lAllen et al.l I2002D . The 

strength of the nucle us destroys any info rmation on the 
color-map. However, ICohen et al.l ()1999l ) found a strong 
Balmer decrement, and a likely V band extinction of 
Ay k, 1.7 mag, u sing K eck spectropolarimetry. They 
and iPrieto et al.l (|1993l ) conclude that the intrinsic 
brightness of the object is likely to put it in the quasar 
class. We would conservatively estimate Ah ~ 0.3 based 
on the se earlier find ings, assuming a diffuse ISM dust 
model ()Mathislll990T ). 

3C 314.1 is an abnormally radio-powerful FR I 
radio source (log 10 Lns = 25.42) in a faint host 
galaxy { Mh = —22.96). It is appar ently not in any 
cluster f Arli ngton- Smith et al.l I1993D . exhibiting an 
under-density of nearby galaxies relative to an average 
elliptical. It shows no major extinction features in 
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elliptical with M H (Host) = -22.99 at z = 0.03. There 
are no signs of disturbance, with ve ry circular (e ~ 0.04 ) 
isophotes on the optical image ([Martel et all 11999ft . 
The color-map shows a smooth gradient in extinction, 
with the highest values up to V = H = 2.7 at the 
cente r, where a small dust lane is visible near the 
core (|de Koff et al.ll2000ft . We estimate A H « 0.1 mag 
for this source i n the central region, assuming a diffuse 
dust model after lMathislll990l . 

3C 402 is an FR II radio source hosted by a smooth, 
elongated elliptical galaxy with M#(Host) = —22.94 at 
z = 0.023. The galaxy is found to be quite red, with 
V — H = 3.0 on average, rising to around 3.5 in the 
central regions. 

3C 445 is an FR II radio source hosted by a 
Mfj(Host) = -22.21 round elliptical at z = 0.056. 
In both the optical and IR images, the morphology is 
dominated by the bright nuclear point source. This 
point source dominates the color-map whi ch thus pro- 
vides no useful information on extinction. ICohen et al.1 
(1999) deduce between about 1.1 and 2.5 mag of V 
band extinction from their spectropolarimetric study. 
This is enough provide up to 0.45 H mag of ext inction 
assuming a diffuse dust model after lMathisl[l990l 

All of the sources discussed above merit closer atten- 
tion to determine their stellar and black hole masses. 
We note that all except two (3 C 61.1, 3C 445)) of these 
low-mass objects are above the iKauffmann et al.l (2003) 
3 x 10 10 M Q cutoff (given our mass-to-light ratio assump- 
tions) but all are significantly less luminous than typi- 
cal quasar host galaxies. Two of the abnormally faint 
FR II host galaxies (3C 61.1, 3C 198) are found to be 
the brightest galaxies in a small cluster or group. The 
remaining six FR II's show anomalously red colors, but 
insufficient evidence of extinction to bring them up to 
L*. The FR I source 3C 314.1 is the only unusually blue 
galaxy. None of the faint objects correspond to the five 
disky (\ow-n) FR II galaxies discussed in Paper II and 
mentioned in Section |4~T1 Those objects were all found 
to have close companions and either to be undergoing 
merger or to be likely recent post-merger systems. 

Other published exceptions to the general rule in- 
clude PKS 0131-36, 0313-192 and PKS 1413+135. 
PKS 0131-36 is a nearby FR II radio source hosted 
by a typical SO galaxy, NGC 612 (see lEmonts et al.l 
120081 and references therein). 0313-192 is a powerful 
(~ 10 24 /if 5 2 W Hz" 1 at 20 cm) FR I ra dio source in 
a highly fla ttened disk-dominated galaxy (|Ledlow et al.l 
[19981 12001 . PKS 1413+135 is a "red quasar" with 
BL Lac object properties in the NIR, but no notice- 
able optical nucleus, located in a disky galaxy with l arge 
amounts of dust obscuration ([Perlman et al.l |2002| and 
references therein). 

6. QUASARS AND RADIO GALAXIES IN CONTEXT 

It remains curious that the low-luminosity and low-n 
sources do not appear in existing samples of RLQ host 
galaxies. This may be simply explained through the 
different methods of selection: in the 3CR, we observe 
all the brightest radio-objects in the northern sky, and 
these may be bright due to any combination of intrinsic 
lumin osity, proximity and beaming. In the R LQ sam- 
ples of lDunlop etHl ([2003ft ; iFlovd et al.l (|2004h . care has 



been taken to remove artificially "radio-loud" sources by 
avoiding sources that are likely to be beamed - objects 
that would appear in the 3CR are drop ped from the RLQ 
samples. However, we also note that iBest et al.l ([2005ft 
find no correlation between radio power and black hole 
mass for radio-loud objects. Thus, it should not necessar- 
ily be expected that beamed sources should have fainter 
host galaxies. 

Much of the drive in quasar host galaxy studies has 
been to confront the technical challenge of observing the 
host galaxies most luminous quasars. In selecting for op- 
tically luminous quasars (radio-loud or radio-quiet) there 
is a dual selection effect for larger, more luminous host 
galaxies: A priori we do not know whether the lumi- 
nosity is host galaxy or AGN-do minated; and if opt ical 
quasars push the Eddington limit (jFlovd et al.ll2004h we 
are likely select for more massive black holes and thus 
(through the well-known black hole to bulge mass corre- 
lation) more massive host galaxy bulges. Radio selection 
does not necessarily select for high black hole mass since 
the radio luminosity is 2 -3 orders of magnitude low er 
than the Eddington limit (jGhiserlini and Celottll2001l) . 

It is possible that samples of RLQs are simply not suf- 
ficiently large to have detected such comparatively rare 
faint objects, of which we only find eight in the present 
3CR sample, or that techniques are insufficiently so- 
phisticated to detect them with current technology (i.e., 
"naked quasars"). Additionally, it is possible that quasar 
host galaxy luminosities can be overestimated du e to 
scattered light from the nucleus (lYoung et alJl2009ft . or 
due to PSF contamination issues (|Kim et al.l l2008). Wc 
feel that this latter effect is unlikely to be a problem 
in the Dunlop, Taylor and Floyd samples, given the de- 
tailed st atistical analys i s and t esting of the method pro - 
vided bv lMcLure et all (fl999l ): lFlovd etaLl (|200l [20081) . 
Finally we note that any strong star formation in the 
quasar host galaxies will have the effect of boosting the 
optical luminosities relative to the simple model assumed 
here. This will in effect bring the bulk luminosities of the 
quasar host galaxies down toward the mean luminosity 
of our sample. 

While there is a noticeable disagreement with the re- 
sults of quasar host galaxy studies, we find that the 
3CR radio galaxies are not unusual with respect to the 
normal elliptical galaxy population. Indeed, only two 
FR II sources - 3C 61.1 and 3C 445 - have galaxies dim 
enough to fall below 3 x 10 *° Mp , ([Kauffmann et al.ll2003t 
iKhochfar andSfil [2006L [2009) assuming the simplis- 
tic co nversion from H band luminosity of IZibetti et al.l 
(2002). These two objects likely have sufficient (~ 
0.5 mag) H band extinction to bring them above the 
mass cutoff without assuming a very different mass-to- 
light ratio. Thus, we have a dilemma: powerful radio 
sources appear to be able to form in even modestly mas- 
sive elliptical galaxies (albeit in low numbers), whereas 
RLQs require a much more luminous host galaxy. More- 
over, studies of RLQs and RQQs have generally found 
a requirement for RLQ host galaxi es to be more lumi- 
nous than their R QQ counterparts (jDunlop et al.l [20031 : 
IFlovd et ail 120041 ) . Thus, selecting for both optical and 
radio luminosity results in more luminous host galaxies 
than both selecting for optically luminous objects and for 
radio-luminous ones. Th ere is no obvious r eason why ei- 
ther the scattered light ([Young et al.ll2009ft or PSF con- 
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tamination (|Kim et al.l 2008) effects should be more sig- 
nificant in RL Qs than in RQQs . We therefore echo the 
conclusion of iBest et al.l (|2005h that radio and optical 
AGNs are separate phenomena. 

7. CONCLUSIONS 
Here, we again summarize the main findings. 

• In terms of ellipticity, the low-z 3CR ar e found 
to be in excellent agreement with iPahrd ()1999f ) 
sample of ellipticals drawn from across a range of 
environments, and with Floyd and Dunlop quasar 
samples. The mergers within the 3CR (identified 
as such in Paper II) fit in well with the genera l 
merger population of iRothberg and Joseph! (|2004[ ) . 

• In terms of Sersic index, there is good agree- 
ment with the Floyd et al. quasar and Urry et al. 
BL Lac object samples, and very poor agreement 
with the merger population, even when considering 
just the 3CR mergers in isolation from the remain- 
der of the 3CR sample. 

• In terms of host galaxy lumino sity, the 3CR 

are g e nerally well matche d to the iBender et al.l 
jlpl); IPahrd (fT999l) and iMobasher et"aTT Tl999) 
elliptical galaxy samples, but exhibit far more faint 
objects than the quasar samples and the BCG's. 

• The Virgo cluster early- type galaxy sample of Fer- 
rarese et al. offers an interesting contrast to the 
3CR, clearly illustrating that the latter is drawn 
from only the most luminous section of a cluster's 
population. 

• The 3CR and its merging subsample have 
a similar luminosity range to the mergers 
of IRothberg and Joseph! |2004), but the merging 
radio galaxies are generally found to be larger, with 
the Rothberg mergers on the low side of the Kor- 
mcndy relation. 

• In terms of radio-optical properties, the FR I's in 
the 3CR unify well with the properties expected of 
BL Lac objects, and FR IPs with RLQs. However, 
a larger spread is seen in both the morphology and 
the host galaxy luminosity of the 3CR sources than 
would be expected from samples of BL Lac objects 
and RLQs studied so far. 

We confirm findin gs fr om earlier work 

(e.g., iLedlow and Owed 119961 : iDunlop et all l2003| ) 
that an elliptical host galaxy is a prerequisite for 
radio-loud AGN activity. However, we identify several 
radio galaxies (~ 20%) that fall below the observed 
RLQ host galaxy luminosity cutoff at ~ L*. These 
objects have a luminosity distribution closer to that of 
the normal elliptical galaxy population, and were missed 
in the NVSS-6dFGS survey due to the K band flux 



limit. The same finding is echoed in the morphological 
comparison - the 3CR host galaxies exhibit a range of 
Sersic index that is entirely consistent with the general 
giant elliptical galaxy population, whereas quasar host 
galaxies show a narrower range of morphology. 

We conclude that the morphological difference is due to 
quasar selection effects and / or contamination of quasar 
host galaxy flux by scattered nuclear flux. It would seem 
highly probable that given accurate coronographic ob- 
servations, or a higher resolution image of the central 
regions of quasars, we would observe the same spread in 
morphologies in the host galaxies of these objects as well. 
However, we do not believe that the luminosity discrep- 
ancy is so easily explained. The fainter FR II host galax- 
ies identified here clearly merit more detailed study and 
form an interesting subset of the 3CR, containing dusty 
sources and brightest apparent members small groups. 

While we now have a reasonable view of the quasars of 
host galaxies, (and a n impressively deep v iew in a small 
number of cases - see lBennert et al1l2008h . several prob- 
lems remain. One is the technical issue of separating any 
scattered nuclear flux from the host galaxy flux. Note 
that this is distinct from the problem of simply separat- 
ing out the host galaxy f r om th e PSF. The problem is 
reviewed by lYoung et alj (|2009T ) who argue that the ef- 
fect is likely to affect the measured luminosities of quasar 
host galaxies significantly. The extent of the problem can 
be easily tested using space-based optical polarimetry of 
a bright quasar. The second issue is one of sample bias. 
Existing quasar host galaxy studies focus on small sam- 
ples with likely selection biases toward bright host galax- 
ies. We need a statistical survey of low-redshift quasar 
host galaxies in at least two bands in order to be able 
to place meaningful constraints on the masses of these 
objects to provide a baseline against which to compare 
higher redshift samples. 
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TABLE 2 

Full sample statistics: D and p from the two-sided K-S test comparing the full 3CR sample with each 

OTHER SAMPLE. 



Property 


U+00 a 


MM01 6 


D+03 c 


F+04 d 


T+96 e 


RJ04-f 


BBF929 


P99 h 


1^+99* 


F+06-' 


G+96 fe 


M H (Rost) 


0.47 


0.34 


0.62 


0.57 


0.44 


0.29 


0.23 


0.15 


0.11 


0.68 


0.81 




0.00 


0.06 


0.00 


0.00 


0.00 


0.01 


0.00 


0.11 


0.82 


0.00 


0.00 


R \/2 


0.33 


0.32 


0.52 
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0.62 


0.30 
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0.26 


0.32 


0.65 


0.60 
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0.00 
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e = l- b/a 


0.36 


0.15 
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0.22 


0.26 
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0.06 




0.20 




0.00 
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0.00 


0.00 
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0.27 
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0.34 
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0.73 




0.00 








0.00 


0.00 


d 












0.42 


0.28 






















0.00 


0.00 











References. 



a lUrry et all {3OT]); " [McLeod e 
^Rothbcrjr and Joseph 2004); S Bender et al. (1992); 



id McLeodl 120011) ; HDunlop et alj liSOOa ); " tt'loyd et"aT1 120041 ); HTaylo 
IPahre] 119991 ); IMobasher et al.| 119991) ; - IFerrarese et al.l 120061) ; ^Ciraham" 



] 119961 ); 
119961 ) 



TABLE 3 

FR I STATISTICS: D AND p OF THE TWO-SIDED K-S TEST COMPARING THE 3CR FR I SAMPLE WITH EACH OTHER 

SAMPLE. 



Property 
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TABLE 4 

FR II STATISTICS: D AND p OF THE TWO-SIDED K-S TEST COMPARING THE 3CR FR II SAMPLE WITH EACH OTHER 

SAMPLE. 
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